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Abstract 

The origin and evolution of isotopes of the lightest elements H 2 , He 3 , Li, 
Be, B in the universe is a key problem in such fields as astrophysics of CR, 
Galactic evolution, non-thermal nucleosynthesis, and cosmological studies. One 
of the major sources of these species is spallation by CR nuclei in the interstellar 
medium. On the other hand, it is the B/C ratio in CR and Be 10 abundance which 
are used to fix the propagation parameters and thus the spallation rate. We study 
the production and Galactic propagation of isotopes of elements Z < 5 using the 
numerical propagation code GALPROP and updated production cross sections. 

1. Introduction 

It has been shown recently [9] that accurate measurements of the antipro- 
ton flux made during the last solar minimum [12] pose a challenge to conventional 
CR propagation models. A solution has been proposed [10] that the observed CR 
may contain a fresh, local, "unprocessed" component at low energies (LE), per- 
haps associated with the Local Bubble (LB). This component reduces the produc- 
tion of B at LE allowing to fit the B/C ratio and antiproton flux simultaneously. 
For more discussion see [11]. 

Isotopes of the lightest elements H 2 , He 3 , Li, Be, B in CR are almost 
all secondary, being produced in spallations of heavier nuclei on interstellar gas 
Galaxy-wide. Their abundances and isotopic composition, therefore, might help 
to distinguish between the propagation models. In this study we calculate abun- 
dances and isotopic composition of H, He, Li, Be, B in two models, a conventional 
reacceleration model [9], and a model with a local component [10]. 

2. The Model 

In our calculations we use the propagation model GALPROP (2D option) 
as described elsewhere [9,14]. The nucleon injection spectrum of the Galactic CR 
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Table 1. Propagation Parameter Sets 





Injection 


Diffusion coeff. at 3GV 


Alfven speed 


Source 


Model 


index, 7 


D , cm 2 s" 1 Index, S 


va, km s _1 


abundances 


A 


1.94/2.42 


6.25 x 10 28 0.33 


36 


LE=HE 


B 


1.69/2.28 


3.30 x 10 28 0.47 


23 


LB^Galactic CR 



is taken as a modified power law in rigidity [6], for the injected particle density. 
The LB spectrum is taken as a power law in rigidity with exponential cut off. 
The Galactic halo size is fixed at 4 kpc. To reduce possible errors due to the 
cross sections we use our own fits to the data on reactions p + He, C, N, O — > 
H 2 ,He 3,4 and p + C,N,0 — > Li,Be,B that produce most of these elements [8]. 
The heliospheric modulation is treated using the force-field approximation; here 
we use $ = 500 MV for all plots. 

In the case of a conventional reacceleration model (model A) we use our 
standard methodologies: the propagation parameters were derived from the fit to 
the B/C ratio, while source abundances were tuned to ACE data [15] and HEAO-3 
data [4] at high energies (HE). No special tuning was done for antiprotons. 

In the case of the model with a local component (model B) we use the 
following procedure. The HE part of the B/C ratio plus antiproton flux mea- 
surements are used to restrict the value of the diffusion coefficient and its energy 
dependence, while the LE part of the B/C ratio is used to fix the reacceleration 
level and define the parameters of the LB component. In this way the model 
provides the best fit to all data at the cost of extra free parameters. 

Model parameters are given in Table 1; model A is similar to models 
DR/DRB as defined in [9], and model B is similar to model DR II as defined in 
[10]. To better match the He spectrum (Fig. 1), a major contributor to H 2 and 
He 3 production, in both models we introduced a break (at 14 GV in A, and at 
10 GV in B) in He 4 (only) power-law injection spectrum. 

3. Results 

Fig. 2 shows the H 2 /He 4 ratio. The measured ratio [2] is a factor of ~1.8 
larger than expected. While there are important hints that this may be caused 
by systematic effects, the analysis is still in progress (D. Vasilas private comm.). 

Fig. 3 shows other ratios. Both models are consistent with data on isotopic 
ratios of Li, Be, and B given the large error bars. He 3 /He 4 ratio is more sensitive. 
Both models agree well with the data, however to match the He 3 /He 4 ratio and 
He spectrum, model B requires a factor of ~2 larger LB abundance of He 4 . The 
HE CR source abundance in both models He/Si w 100, while the LB abundance 
in model B is 220 (cf. solar system value 2400). Taking into account that He and 
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Kinetic energy, GeV/amu Kinetic energy, GeV/nucleon 



Fig. 1. He spectrum. Calculations: ••• 
- model A, B. Solid line - mod- 

ulated (A k. B), thin lines - interstel- 
lar (LIS). Data: for references see [9]. 



Fig. 2. H 2 /He 4 ratio. Calculations: • • • 

- model A, model B. Bold lines 

- modulated, thin lines - local inter- 
stellar ratio (LIR). Data: IMAX92 [2]. 



Si are abundant elements a factor of 2 difference is significant. Contributions of 
C 12 and O 16 to He 3 and He 4 production appear to be non-negligible. 

The smaller proportion of He compared to heavier elements in Galactic 
CR sources supports the idea that HE CR are accelerated in shocks from fresh 
SN ejecta. The higher proportion of He compared to heavier elements in model 
B at LE implies that the material was diluted before acceleration. This is in line 
with the view described in [10] that the LB component may be accelerated by an 
ensemble of weak shock waves out of the interstellar medium. If the model E is 
correct it may indicate that the LE part of the He and proton Galactic spectra 
is flatter than thought with corresponding consequences for the diffuse Galactic 
7-ray emission and Galactic chemical evolution. 

This work was supported in part by a NASA Astrophysics Theory Program 
grant and by the US Department of Energy. 
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Fig. 3. Isotopic ratios as calculated in two propagation models. Lines are coded as in 
Fig. 2. Data: IMAX92 [13], ISOMAX98 [5] (with statistical errors only), Voyager 
[7], Ulysses [1], ACE [3]. 
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